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Abstract 

In this paper we show how a change of box dimension of the orbits of 
two-dimensional discrete dynamical systems is connected to their bifur- 
cations in a nonhyperbolic fixed point. This connection is already shown 
in the case of one-dimensional discrete dynamical systems (see [12j.[8]l. 
Namely, at the bifurcation point the box dimension changes from zero to 
a certain positive value which is connected with the type of bifurcation. 
First, we study a two-dimensional discrete dynamical system with only 
one multiplier on the unit circle, and get the result for the box dimension 
of the orbit on the center manifold. Then we consider the planar discrete 
system undergoing a Neimark-Sacker bifurcation. It is shown that box di- 
mension depends on the order of the nondegeneracy at the nonhyperbolic 
fixed point and on the angle-displacement map. As it was expected, we 
prove that the box dimension is different in rational and irrational case. 
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1 Introduction 



In recent years, it was shown that fractal analysis can be applied to the solu- 
tions of differential equations and dynamical systems (see [17], [20], [21], [30]V 
In the case of dynamical systems, fractal analysis consists of studying the box 
dimension and Minkowski content of trajectories or orbits. Several articles with 
fractal analysis of bifurcations of dynamical systems (see [5], [21], [IH], [29) ) 
showed that there is a direct connection between the change in box dimension 
of trajectories of dynamical systems and the bifurcation of that system. Around 
the hyperbolic singularities the box dimension is trivial (0), while around the 
nonhyperbolic singularities the box dimension is positive and connected to the 
appropriate bifurcation. In fact, box dimension shows the multiplicity of singu- 
larity or limit cycle in the case of weak focus or multiple limit cycle (see [2H])- 
Beside this, in the article [18], multiplicity was also connected to the growth of e- 
neighbourhood of the orbit near the homoclinic loop. At the beginning, this phe- 
nomenon was studied only in the continuous dynamical systems, but the same 
occur in the discrete systems. The article [8^ showed that the above result for 
box dimension can also be proven for generic saddle-node and period-doubling 
bifurcations of one-dimensional discrete dynamical systems. In the paper [12j 
the result from the previous paper were generalized to the class of finitely non- 
degenerate maps in R and applied it to one and two-parameter bifurcations 
with generalized sufficient conditions. So in the one-dimensional systems, the 
box dimension at the bifurcation point shows the maximum number of possible 
fixed point that can bifurcate in a given family of systems. Of course, now we 
are interested in this connection for other bifurcations. 

Hence, it is natural to continue with the fractal analysis of two-dimensional 
discrete dynamical systems. For the result about the trivial box dimension of 
the orbit near the hyperbolic fixed point in E", see the article [13]. In this 
article we will show the results about the box dimension of the orbit around 
a nonhyperbolic fixed point in M^. We look at the following bifurcations: one 
multiplier on the unit circle and two complex conjugated multipliers with the 
nonresonant condition. For the bifurcation of nonhyperbolic fixed point with 
only one multiplier on the unit circle, we get similar result for the orbit on 
the center manifold, which correspond to the bifurcation of appropriate one- 
dimensional system. Regarding the Neimark-S acker bifurcation (see [19], [24], 
|15|). the situation is more complicated. As in the bifurcation analysis of this 
bifurcation, the main problem is the difference of the dynamics between the 
irrational and rational case. We will see that the box dimensions for this cases 
are also different. Since Neimark-Sacker bifurcation is, in fact, Hopf bifurcation 
for maps, we will see how this result is related to the result for Hopf bifurcation 
showed in [28] . 

The main motivation for studying the box dimension of dynamical systems 
is related to the problem of multiplicity or cyclicity of fixed point, singularity 
or limit cycle. We are interested to explore how box dimension behaves during 
the bifurcations in of discrete and continuous systems, its connection to the 
cyclicity problem and how it can be found helpful in bifurcation analysis of some 
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complicated bifurcations in higher dimension. Moreover, this connection could 
be applied as numerical tool in bifurcation analysis since there exists effective 
algorithms for calculating box dimension. 

Now, we will introduce the notions of box dimension and Minkowski con- 
tent, as basic terms which are used in our fractal analysis. First, we recall that 
the fractal dimension such as box dimension (also known as Minkowski dimen- 
sion, Minkowski-Bouligand dimension, capacity dimension, limit capacity), can 
be used to analyze various objects such as various sets, graphs of a function, 
attractors, trajectories, etc. Fractal analysis consists of getting the results of 
the fractal dimension and putting them in the context according to other prop- 
erties of the studied object. In the case of dynamical systems, we connect the 
box dimension with the bifurcation, its type and the number of bifurcating ob- 
jects. In our analysis we use the box dimension because, for the orbit of discrete 
dynamical system, the Hausdorff dimension fails to show anything. Namely, 
because of its property of countable stability, the Hausdorff dimension does not 
'see' the countable sets at all. On the other hand, the box dimension is only 
finitely stable so it can 'see' the countable sets and the difference between the 
orbits around the hyperbolic and nonhyperbolic fixed point. 

Now we recall the notions of box dimension and Minkowski content. For 
further details see e.g. [H] , Hi] , [17] . Let C be bounded. The e- 
neighbourhood of A is defined by = {y ^ : A) < e}. 
Let s > 0. The lower and upper s- dimensional Minkowski contents of A are 
defined by 

:=liminf limsup 

Then the lower and upper box dimension are defined by 

dimpA = inf{s > : MliA) = 0}, dhRgA = mi{s > : M*'{A) = 0}. 

If diui pA — dims A we denote it by dims A. If there exists d > such that 
< A4'^{A) < A4*'^{A) < oo, then we say that set A is Minkowski nondegener- 
ate. Clearly, then d ~ dims A. li\A^ \ ~ e'^ for e small, then A is Minkowski non- 
degenerate set and diuiB A = N - s. If Ml{A) = A^"(A) = M'''{A) e (0,oo) 
for some d > 0, then A is said to be Minkowski measurable. Clearly, then 
d = dims A. 

Let A C M" be a disjoint bounded set and F : R" ^ R" is a Lipschitz map. 
Then it holds 

dims F{A) < dims A. 

We say that F : ^ fl' , where il, fl' are open sets, is a bilipschitz map if there 
exist positive constants A and B such that 

A\\x~y\\< \\F{x)-F{y)\\<B\\x~y\\, 

for every x,y €z If -F is a bilipschitz mapping, than 

dims A = dims F{A). 
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Box dimension has also the property of finite stability: Let Ai, i = I, . . . , k 
be disjointed bounded sets. Then 



If dims Ai = d for every i, then dims y_j Ai = d. 

In the paper the following definitions are also used. We say that any two 
sequences (a„)„>i and (&«)«>! of positive real numbers are comparable and 
write a„ ~ 6„ as n — > oo if j4 < a„/6„ < B for some A,B>Q and n sufficiently 
big. Analogously, two positive functions f,g : (0,r) — > M are comparable and 
we write f{x) ~ g{x) as x — !■ if f{x)/g{x) e [A, B] for x small enough. 

We will study a two-dimensional discrete dynamical system 



generated by a function F : i— >■ . The orhit of a system is a sequence 
(x„)„>i such that x„+i = F(x„) for some xi S K^. Let cco = be a fixed point 
(F(xo) = xq) of that system and let A be a Jacobi matrix DF{xq) at xo- The 
eigenvalues Ai, . . . , A„ of the matrix A are called the multipliers of fixed point. 
We denote by no the number of multipliers on the unit circle, by n_ the number 
of multipliers inside the unit circle and by the number of multipliers which 
lies outside the unit circle. The fixed point is hyperbolic if no — 0, that is, 
there is no multipliers on the unit circle. Hyperbolic point is called a hyperbolic 
saddle if ^ 0. The fixed point is nonhyperbolic if uq ^ 0. We will also 

need the following definition (see |12)). 

Definition 1 Let F : {xq — r, + r) — > M, r > 0, be a map of class , and 
Xo is a fixed point of F such that F'{xq) ^0. // there is a k > 3 such that 
F"{xq) = . . . = F^''~^^{xo) ~ and F^^^{xq) ^ 0, then we say that the map 
F is a k-nondegenerate map in xq. Specially, if F"{xo) ^ 0, then we say 
that F is 2-nondegenerate map in xq. The number k is called the order of 
nondegeneracy of map F in xq. 

In this paper the main object of our study is a box dimension of the orbit 
around the nonhyperbolic fixed point of discrete planar dynamical system. In 
Section 2, we prove the result for the two-dimensional systems with only one 
multiplier on the unit circle by using the center manifold theory. In Section 3, we 
analyze the Neimark-S acker bifurcation. We will see that there is a difference in 
box dimension in rational and irrational case. The rational case will be proven 
by direct calculation, while the irrational case is showed using the analogous 
result for the Hopf bifurcation from the article [28] . 

2 Fractal analysis of bifurcations in M? 

We will begin the study of bifurcation of discrete dynamical systems in with 
the bifurcations of two dimensional maps with only one multiplier on the unit 
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(1) 



x„+i = F(x„), xi e 
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circle. Without loss of generality, the system can be written in the form 

Xn+1 = MXn+ liXn.Vn) (2) 
Vn+l = >^2yn+ 9{Xn,yn) 

where |Ai| = 1 i IA2I < 1, and /, g are of class C on some neighbourhood 
around the origin such that /(O, 0) = 0, £1/(0, 0) = 0, 5(0, 0) = 0, £15(0, 0) = 0. 
It is clear that (0, 0) is a nonhyperbolic fixed point of the system. All other cases 
for multipliers can be proven analogously. Then by Center Manifold Theorem 
we know that there exists C""-center manifold y = h{x), and the restriction of 
the system is 

Xn+1 = \lXn+ f{Xn,h{Xn)) (3) 
yn+l = h{Xn+l). 

Now we would like to determine the connection between the box dimension 
of the restriction and its projection on the x-axis which is 

Xn+1 = \iXn + f{Xn,KXn)) (4) 
yn+l = 0. 

We know that the projection on the a;-axis is a Lipschitz map, but in order to 
get the equality of the dimensions we need the bilipschitz mapping. At this 
point, we will use the next lemma from the article [29] . 

Lemma 1 (see ]2ff]) Let g : R^^^ R be a Lipchitz map, N > 2, and we 
define F -.R^ ->-R^ with F{x, z) (x, z + g[x)), where x € R^-^ and z G R. 
Then F is a bilipschitz map and holds the measure, that is, for every measurable 
set E CR^ of limited measure it holds = \E\. Furthermore, for every 

limited set A C R"'^ we have 

dims F (A) = dims A i dim pFfA) = dim pA. 

Set A is a nondegenerate if and only if F{A) is a nondegenerate. 

We apply this lemma on the plane {N = 2) and can see that if g : R ^ R is 
a Lipschitz map, then the map F : R^ — R^ defined by 

F{x,z) = {x,z + g{x)) (5) 

for x,z € R is a bilipscitz map. 

Now by using the above lemma we can prove the theorem about the box 
dimension of a discrete system on the center manifold. 

Theorem 1 Let the restriction of the system ^ with Ai = ±1 on the center 
manifold y = h{x) is given by 

x I— >■ G{x) = Xix + f{x, h{x)). 
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Let S{xi, yi) = {x„, Un) be an orbit of the systems on the center manifold in the 
form 

Xn+l = XlXn + f{Xn,h{Xn)) (6) 
Vn+l = HXn+i) 

with initial point {xi,yi) near (0,0). // the map G is a k-nondegenerate in 
X = 0, then there exists r > such that for \xi \ < r we have 

dims S{xi,yi) = 1 ^ 

Proof. First we define the set A — {{x,y) : x € A{xi),y = 0}, where 
A(xi) = {xn)ni£N IS a one-dimensional discrete dynamical system generated by 
Xn+i = G(xn) — Aix„ + f{xn,h{xn)) and xi € (0,r). Then we act with the 
map ([5]) on the set A, i.e. with F{x, z) — {x, z + h{x)) and get 

F{A) = F{x, 0) = (x, + h{x)) = {x, h{x)). 

So the image of the set A under map F is associated system ([6]) on the center 
manifold y = h{x). In other words, the map F associate the projection of a 
orbit on a;-axis with appropriate orbit on the center manifold. Since the map h 
is of class C on some neighbourhood small enough < 6 and h'(0) = 0, then 
h' is limited on that neighbourhood, and we have 

\h{x2) - h{xi)\ = \h'{x*)\ \X2 -Xi\< M \X2 - Xi\ 

for some x* £ {xi,X2). Therefore h is a Lipschitz map for |a;| < 6. Now it 
follows from Lemma 1 that is a bilipschitz map, and 

dims F{A) = dims A. 

Since G is a fc-nondegenerate map, from [T^] (Theorem 2.2), it follows that 
there exists r > such that for the sequence A{xi) — (x„)„gN defined by 
Xn+\ = G{xn), \xi \ < r we have 

dims A{xi) = 1 ~ 

We see now that dims A = dims A{xi) — 1 — ^. Notice that S{xi,yi) — F{A), 
so it follows that for jxij < r is 

dims S(xi, yi) = dims F(A) = dims A ~ \ — — . ■ 

k 

Remark 1. The previous result can be easily generalized to the dynamical 
systems in R" with only one multiplier on the unit circle. 

Next we will apply the previous result to the systems with Ai = lorAi=— 1. 
By short calculation we get the following forms of the restrictions and center 
manifolds. 



5 



The center manifold of the system ^ with Ai = ±1 is given by 
V = V{u) ^ ^uju^ + 0{u^), 
and the restriction of the system ([2]) on that center manifold is 

where ^, cr, a, 6 G K are specified by 

^_ ^^(0,0) g^/(o,o) ^_ 5^g(o,o) ^_ 5^(0,0) 

9a;2 ' dx^ ' cJa::^ ' dxdy 

Notice that in the cases of saddle-node, pitchfork, transcritical and period 
doubling bifurcations the box dimension can be established from ([7]). For details 
see [12], [15], [23]. So, in the case of saddle-node and transcritical bifurcation 
(7 7^ 0, so the box dimension of a orbit of the system on the center manifold for 
Xi small enough is 

dims S{xi,yi) = -. 

For pitchfork and period doubling bifurcation is a = 6 = and 5 7^ so we have 

2 

dims S{xi,yi) = -. 
Now we can look at some examples. 

Example 1. We consider the family of discrete planar systems of a form 

X i-> /i + X + (8) 

y ^ y/2 

with parameter fi gM.. Notice that a = so wi = and the centre manifold is 
V = 0. From a = 2 and (5 = 0, it follows that the restriction of the system for 
/i = on the center manifold v = is of a form 

As we can see at Figure 1, the restriction of the system on the center manifold 
exhibits saddle-node bifurcation, and for fi — the fixed point (0, 0) is nonhy- 
perbolic with the box dimension dimes' = ^. After the bifurcation we have 
two fixed points on the x-axis. Also notice that the box dimension of orbits of 
one-dimensional system y 1-^ y/2 on y-axis near = is 0, since ?/o = is a 
hyperbolic fixed point for that map. So we can say that we have hyperbolic and 
nonhyperbolic direction. 
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Figure la fi <0 



Figure lb jj, = 



f " — ji ■ 

■ . i . ' ! 

Figure Ic > 



Example 2. We consider the family 

X H> (/u - + (9) 
2/ y/2 

with the parameter fj, gR. Notice that a = so W2 = and the center manifold 
is w = 0. Also, notice that cr = so the restriction of the system for = on 
the center manifold w = is of a form 

u —u + u^. 

As we can see at Figure 2, the restriction of the system on the center manifold 
exhibits the period doubling bifurcation, so for = the fixed point (0,0) is 
nonhyperbolic with box dimension dims 5' = |. 
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Figure 2a fi < 

T 

• ■ - - - - , ~ 

r - - ^» ^ . 



Figure 2b fi ^ Figure 2c fi > 

Considering the bifurcations with only one multipher on the unit circle, we 
see that in this cases only more fixed or periodic points can bifurcate from one 
nonhyperbolic fixed point. In order to get the objects of bigger dimension from 
the bifurcation (e.g. invariant circle), we must have at least two multipliers on 
the unit circle. 

The above-mentioned procedure for getting the box dimension of the orbit 
on the center manifold can be analogously applied to the orbits on the stable 
and unstable manifolds (see [H]). This result can also be applied to appropriate 
bifurcations of continuous dynamical systems (see |13)). 



3 Fractal Analysis of Neimark-Sacker bifurca- 
tion 

In this section we study the bifurcation at the nonhyperbolic fixed point with two 
complex conjugate eigenvalues on the unit circle when invariant curve bifurcates 
called Neimark-Sacker bifurcation. For details, see [TS], |19 ) .|24 ) . 

From |15| . we know that the polar normal form of one-parameter family of 
maps for Neimark-Sacker bifurcation is 

r r + d^r + ar^ + 0{^i'^r, ^r^,r*) ^ F{r) 

if ^ if + ea + QiH + br^ +0{n'^,fir\r^)^G{^,r) (10) 
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where S R, and coefficients d, a, 0o,6i,6 € M satisfy tlie bifurcation condi- 
tions: a ^ 0, d ^ and e'''^° ^ 1 for fc = 1,2,3,4. Since tlie map F doesn't 
depend on tlie angle (p, studying the existence of the invariant curve is simple. 
On the other hand, the map G depends on r, which complicates the study of 
the orbit structure on the invariant curve, and consequently, the calculation 
of box dimension. All the previously studied bifurcations (see [5], [5S], [H]) 
showed that the box dimension of a orbit around the nonhyperbolic fixed point 
is connected with the box dimension of the invariant set which emerge at the 
bifurcation point. For instance, in the bifurcation of one-dimensional discrete 
dynamical systems when only fixed point can bifurcate, the box dimension is 
between and 1. Furthermore, the box dimension in Hopf bifurcation when we 
have limit cycles is greater then 1. 

Moreover, in all other bifurcations of discrete and continuous systems the 
invariant sets which emerge consist of only one orbit. For example, the limit 
cycle which bifurcates from the Hopf bifurcation has box dimension 1, and 
consist of only one trajectory. In the case of Neimark-S acker bifurcation, the 
originate invariant curve consists of many different orbits, and their structures 
depend on the map G{(p, ro) from pO)) . with r = ro the radius of invariant curve. 

We know that the orbit structure on the invariant curve is depended only 
on the map G{ip,r). Therefore, there is a case G{f,ra) — ip + 0o, 6o = ^ 
— )•. Then (p, q')-cycles emerge by the bifurcation and the box dimension of 
every orbit on the invariant curve is (property of finite stability). In the 
case G((/3, rg) = ip + Oq, 6o = 27r/3, /3 irrational, on the invariant curve every 
orbit is dense and its box dimension is 1. It is known that in general case with 
G{(p, ro) = 1^ -t- 6o -I- Oi/Lt -|- brQ + 0(/^^, /xrQ, Tq) — the orbit structure changes 
by changing the parameter between rational and irrational rotation numbers 
(Arnold tongues, see [l5].[T]). 

In the article [13], it was showed that in the neighbourhood of the (un)stable 
hyperbolic fixed point in K", the box dimension of every orbit is 0. Now we 
want to establish what happens with the box dimension of an orbit around the 
nonhyperbolic fixed point. So there is no need to study whole family of one- 
parameter maps but only the system for the bifurcation value ^ = because 
then the point xg = (0, 0) is a nonhyperbolic fixed point, and the box dimension 
is positive. Our goal is to determine the value of that box dimension. 

Hence, we look at the polar normal form for Neimark-Sacker bifrucation at 
the bifurcation value fi = 

r ^ r + ar^ +0{r'^) = ,f{r) 

ip ^ ip + Qo + br^ + Oir^) = g{ip,r) (11) 

where a 7^ (nondegeneracy condition) and e*® 7^ \/T, for n ~ 1,2,3,4 (non- 
resonant condition). We consider the orbits of the system (fTTj) around the 
nonhyperbolic fixed point xq = (0,0). Let a < 0, and define 

F = {{rk,(pk) ■■ Tk = f{ru-i),<pk = g(rfc_i,(pfc_i), fc e N, (ro,<^o) given}. (12) 

When we draw only one orbit F, we observe that the points are spirally going to 
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the origin, and that is why the orbit F is cahed a discrete spiral. To be precise, 
for a < the discrete spiral is spirally going to the origin, but for a > the 
orbit is spirally going away from the origin. Furthermore, if the rotation angle 
is positive, then the spiral has a positive direction, and otherwise the direction 
is negative. See Figure 3. 




Figure 3a f{r) = r — r^, Figure 3b f{r) = r — r^, 

g{r, If) = (p - ^ -r'^ g{r, <f) = (p + l + r'^ 

3.1 Bounds for box dimension 

We know that for every bounded set in the plane we have 

< dimpF < dmiijF < 2. 

Now we recall that the generalized polar normal form for the Neimark-Sacker 
bifurcation proven in the article |^ (Theorem 5) is 



p 

r I— >■ r — 

m— 1 



^a„r2"+^+0(r2^'+2) = /(r) 

r?j.= l 

P 

^ ^ ip + eo+J2b2rnr^"' + Oir^P+')=gi^,r). (13) 

m— 1 

and the corresponding nondegeneracy condition: 

03 = . . . = a2k-i = 0, a2k+i ^ 0. (14) 
It is clear that a = 2fc + 1 is the order of nondegeneracy of the map /. 
Theorem 2 Bounds for box dimension 

Let T be a discrete spiral defined by il2\) where f and g are defined by the system 
P3p . /// is an a-nondegenerate in xq — 0, then there exist rg > small enough 
such that 

1 - - < dimpF(rn, <^o) < dimBF(ro, <^o) < 2(1 - —). 
a a 
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Proof. 

If / is an a-nondegenerate map in a point xq, it means that 

f'ixo) = ... = f'^-'Hxa) - 0, /("'K) ^ 0. 
In the case when f = f{r) in the system p^ . then 

03 = . • ■ = a2fe-i = 0, a2k+i ^ 
for some fee 1, 2, . . . ,p, so then a = 2k + 1. 

Let A be a sequence of concentric circles of radius rk+i — f{rk), k £ N with the 
centre in the origin. Notice that F C A, so in the case of nonhyperbohc fixed 
point 

dinisF < dims A = 2dimB{rk}keT>! = 2(1 ). 

a 

Of course, we have used the fact that the box dimension of the orbit of one- 
dimensional sequence {r^} generated by an a-nondegenerate map is 1 — ^ (see 
|12j). In this way we got the upper bound for the box dimension of F which is, 
in fact, the best upper bound. 

Now we are interested in the lower bound. Let : M.'^ ^ be a radial 
projection on the x-axis defined by Px{r,ip) = (7^,0). This map is a Lipschitz 
map, so 

dim^F > dims PxiT) = dim_B{rfc}fegN- 

It means that we get the following estimation which holds around the non- 
hyperbolic fixed point 

dimB{rfc}fcgN < dim ^F < dinisF < 2 dim_B{rfc}fcgN, 

that is, 

1 - - < dimpF < dWsF < 2(1 - -). ■ 
a a 

So, in the case of classic Neimark-S acker bifurcation for a = 3 we have 

I <dim»F<di^RF< ^. (15) 

In the next section we will shown that for the fixed and rational rotation 
angle the lower bound is achieved, while in the irrational case, the upper bound 
is achieved. 

3.2 Fixed and rational angle displacement map 

The polar normal form of Neimark-S acker bifurcation at the bifurcation value 
fi = with fixed rotation angle (g doesn't depend on r) 

r ^ r + ar°'+0{r°'+^)^f{r) 

ip (p + Qo^ g{ip) (16) 



11 



where a is odd. If the rotation angle is of a form 



eo 



2'Kp 



(17) 



for some p,q € Z, N{p, q) = 1, we say that the rotation angle is rational. That 
is the special case of Neimark-Sacker bifurcation, when by bifurcation emerges 
the invariant circle with the (p, 5)-periodic cycles on it, that is, every orbit on 
the invariant circle consist of finite number of points, so its box dimension is 0. 
Due to the fact that there exists q such that 

gGo = 27rp, 

the orbit T of the system (fT6)) will be the union of q rays from the origin with 
the one-dimensional sequences of points which tend to origin. On Figure 4 we 
can see the phase portraits of the system 



/(r) = (l + /^)r- 
gir, If) 



^+6 



(18) 



depending on /i. Notice that g = 12, so on every picture the points are on the 
rays, but around the nonhyperbolic fixed point they are evident (Figure 4b). 
On Figure 4c we see that the orbit on the invariant circle are (12)-cyclcs. 




Figure 40, /i < Figure 4b n = Figure 4c fi > 

For the box dimension result we will need following lemma. 
Lemma 2 Let be a system 



r ^ r + ar"+0{r°'+^)=f{r) 

with a > 3. Then the k. iteration of that system is of a form 
f\r) = r + kar°'+0{r°'+^) 



(19) 



(20) 



Remark 2. This technical lemma can be easily proven by using the math- 
ematical induction. 
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Theorem 3 (fixed and rational rotation angle ) 

Let r(ro,<po) be a discrete spiral defined by i fl2]) where f and g are given by a 
system ((TE)) . so let the fixed rotation angle be in a form jTT]). Then for ro > 
small enough we have 

diniBr(ro,^o) = 1 ■ 

a 

Moreover, T{ro,(pQ) is Minkowski nondegenerate. 
Proof. 

Because of the fixed and rational rotation angle, for the discrete spiral with the 
initial point (ro,ipo), it holds 

r(ro,¥>o) = UfroV,(ro,(p.) (21) 

where 

Uro,Vi} = {r(ro),r+'(ro),/''+'(ro),...} 

= {fHr„)J\f^ro))j'^f\ro)),...} (22) 

is a one-dimensional sequence on the ray ip ^ ipi = g^{ip()). We see that the 
discrete spiral Ti{ro, ipi) is in fact generated by the map with the initial point 
f^{ro), i = 0, . . . ,q — 1. From Lemma 2, we have /'(r) = r + qar" + 0{r"^^). 
Now from |12^ (Theorem 2.2), it follows that for rg > small enough 

dimBri(ro,<^.i) = 1 - -. 

a 

Now by using the finite stability property (see Section 1) we get 

dims r(ro, cpa) = dims ri(ro, cpa) = 1 - - (23) 

a 

since all the box dimension on the rays are equal for tq small enough. 
Regarding the Minkowski content, it can be easily seen that 

9-1 

|ri,.|<|r,|<^|r,,,|, 

so 

9-1 

MtiTi) < MtiT) < M*'^iT) < ^X*''(r,) 

1=1 

where d = dimsF = dimsF,;. From |12j (Theorem 2.2) we know that F.^ are 
Minkowski nondegenerate sets, so we get that 

0<Mt{T) <M*'^(r) <oo. ■ 
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3.3 General rational case 



In this section we will show the method for getting the box dimension in general 
rational case. In the previous section we got the box dimension result by using 
the finite stability property and some previously known results. In this section, 
we will show the method of direct calculation by using the definition of box 
dimension. First step is finding the good estimation of iFgl, where T — T{ro, ipo) 
is a discrete spiral and is its e-neighbourhood. 

At the beginning, we will study the calculation method with the fixed rota- 
tion number. Hence, we would like to estimate the area of the £-neighbourhood 
of one orbit of the system (IT51) with a < 0, fixed and rational rotation angle 
80 and for fixed e > 0. Without loss of generality we may assume that a = 3. 
For odd a > 3 the calculus is analogous. In order to do that, we have to look 
in detail what are all the possible overlaps in the mentioned e-neighbourhood. 
Notice that the following analysis is valid in the general situation as well. 

In the case of discrete orbit F, the e-neighbourhood of F, we denote it by 
Fe, is a union of circles around points Ak{rk,ipk) with radius e. The main 
problem is that after some k the circles begin to overlap and then it is difficult 
to calculate its overall surface. We observe two neighbouring points Akirkj^Pk) 
and Ak+i{rk+i,ipk+i), and one point in the next level Ak+qa{rk+qo,Vk+qo) where 
go is the lowest positive number such that go^o > 27r. See Figure 5. We will 
see that in the case of rational rotation number with p = 1 is qoQo = 27r. For 
p > 1 the number go will be such that qoOq = 2pTT. Therefore, we have three 
significant distances whose behavior we must analyze in order to determine the 
way and the order of overlaps. 




Figure 5 



We denote with 

yk ^ d{Ak,Ak+i) = Jrl+rl^-^ - 2rfcrfe+i cos(eo) (24) 



Zk = d[Ak,Ak+q,,) = \Jrl + rl^q^ - 2rkrk+q„ cos(goeo) (25) 
Wk = d{Ak+i,Ak+q„) = \Jrl^i + rl^q„ - 2rk+irk+qo cos((go - l)eo)(26) 
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Hence, we have three possibihties 

Vk < 2e (27) 
zk < 2e (28) 
Wk < 2e (29) 

First, for every overlapping we must show that the sequences {yk}, {zk} and 
{uJk} are decreasing to ensure that m{£) is the lowest natural number for which 
overlapping occurs. So we study the behaviour of the sequences {i/k}, {zk} and 
{wk} as A; — >■ oo. Since rk+q„ appears in the expressions for Zk and Wk, we need 
the expression for rk+qg — f'^°{rk) showed in Lemma [2J 

It is easy to prove the lemma about the behaviour oi {yk}, {zk} and {wk} 
for fixed and rational 9o- 

Lemma 3 Let be a system with a < 0. Then the sequences {yk}, {zk} i 
{wk} defined by 121]); and are decreasing for k big enough, and 

Zk ^ (jo|a|fc~^ (30) 

wkc^yk V2(l-cos(eo))fc-^. (31) 

Moreover, there exists Kq such that Zk < yk for all k > Kq. 

Proof. 

First, we show that {yk} is decreasing (yk+i < yk)'- 
By squaring yk we get 

rl+i + rl+2 - 2rfc+irfe+2 cos(eo) < rl + rl^^ ~ 2rkrk+i cos(eo) 

and 

Tk + rk+2 > 2rfe+i cos(eo). 
In the proof of Theorem 2.2, |12) it was shown that one-dimensional sequence 
of a form rk+i — rk + ar'^ + ©(r^^^), for a < has a decreasing sequence of 
differences, that is, (r^ — rk+i) is decreasing on some small neighbourhood of 
the fixed point (for k big enough). Since the box dimension is also calculated 
on some small neighbourhood of the fixed point, it means that for k big enough 
it holds 

rk+i - rk+2 < rk - rk+i 

and follows 

rk + rk+2 > 2rfc+i > 2rk+i cos(eo). 

So {yk} is decreasing for k big enough. Completely analogously can be proven 
that {zk} and {wk} are decreasing as well. 

By including rk+i = rk + ar^ + 0(r^+i) and rk+q, = + qoar^ + 0(r^+i) 



in the expressions (PH) . (1^ and (1^ we get 

yl = 2(l-cos(eo))r^+0(rr') (32) 

zl = 2(l-cos(goeo))r^+0(r,"+i) (33) 

wl - 2(l-cos((go-l)eo))r2+0(r^+i) (34) 
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Using the known behaviour of {rk} as fc — > cxi (Theorem 2.2, [T2]) 
rfc ~ fc~^ i.e., 3A,B > Ak~^ <rk < Bk^^ 

we obtain 

Vk ^ V2(l-cos(eo))fc-^ (35) 



Zk ~ V2(l-cos(goeo))fc-— (36) 
wk ^ V2(l-cos((go-l)eo))fc-^. (37) 

The condition cos(0o) 7^ 1 is satisfied because of the nonresonant condition 
for the Neimark-Sacker bifurcation. This type of behaviours wih be vahd also 
for 80 = 27r/3, /? irrational Of course, if the rotation angle is rational, then 
QoQo — 27rp. Then it follows 

Zk ^ qo\a\k'^ (38) 
wko^yk V2(l-cos(eo)fc^^. (39) 

Regarding the relation between Zk and yk , some minor problem can occur since 
both depend on qq. That is why we will closely examine their relationship. Let 
A^^Bz > be constants such that 

A^qa \a\ k^^ < < B^Qq \a\ fc^^, 

while Ay, By > are such that 

Ay^2il - cos{eo))k-^ < Vk < By^2{l-cos{eo))k-^. 

We want to have 

Zk < B^qo \a\ fc"^ < Ay ^2(1 - cos(eo)fc~^ < yk, 

that is, 

^ ^ Bzqo \a\ 



Ayy/2{l-coseo) 

In order to obtain ]^_gQs(g)^-) we use Qq = and the Taylor expansion for 
cosine, and get the existence of constants Ci, C2 > such that 



27rVCi - 1 -cos Go " 2C27rV 
It means that k need to satisfy 



k > 



^-K-pyfClAy ' 



Now we have that Zk < yk for k > Kq = 2TTpy/c A 
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Remark 3. For a > 0, we need to observe the inverse map for /, i.e. the 

system 

r ^ r - ar" + 0{r"+^) 

ip 1-^ <y3 - 6o (40) 
and the result from lemma also holds. 

Now we study the case 8o — 27r/3, (3 E Q, but the rotation angle isn't fixed 
but has also the terms of higher order. So we observe the system 

(f ^ ifi + + br°'-^ + 0{r°') (41) 

where 6o = p,q £Z , N{p, q) — I, a < and a odd. 
The overlapping sequences are given by 

yk = d(Ak,Ak+i) = \Jrl +rl^^ - 2rfcrfe+i 003(60) (42) 



Zk = d{Ak,Ak+qa)= Jrl^rl ~-2rkrk+q„co8{ipk+qo-Vk) (43) 



Wk = d(Ak+i,Ak+qo) = Jrl_^i+rl -2rk+irk+qocos{<pk+qo - 'Pk+i). {4:A) 



Lemma 4 Let be a system with a — 3. Then the sequences {yk}- {^k} 

and {wk} defined by i(42|) . (j43p and (1441) are decreasing for k big enough, and 

Zk ~ qoVa^+lPk^i (45) 
wk^yk ^ v/2(l - cos(eo))fc"i (46) 
Moreover, there exists Kq > such that Zk < yk for all k > Kq. 
Proof. 

First we will show that {zk} is decreasing for k large enough. Hence, we have 

4 = rl + rl^g^ ~ 2rkrk+qo cos{ipk+qo - fk), 

and we want to show that Zk+i < Zk for large k. It is easy to show that 
ipk+qo+i - ^k+1 < fk+qo - Vk and it follows that 

- cos(<^fc+,o+i - ipk+i) < - cos{(pk+qQ - Vk)- (47) 

Now recall that the sequence of differences for {rk} is also decreasing so 

rk - rk+i > rk+qo - rk+qo+i- 
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By squaring and reducing we get 

rl - rl+i + rl^g^ - rl^qo+i > '^rkrk+qo - '2rk+irk+q„+i > 

> {2rkrk+qo - 2rk+irk+q„+i)cos{(pk+qo - (fk), 

and by using (|47)) we get 

rl+rl+q„-2rkrk+qo coa{(pk+qo-'Pk) > rl^l+rl^q^_^_^-2rk+lrk+qo+lCOs{(pk+qo+l~^k+l)■ 

Analogously for {wk} and {yk}- 

Now we observe the behaviour of Zk ■ In the expression 

4 = rl+ rl^q,, - 2rkrk+qo cos{ipk+q,, - (fik) 

we put 

/9o(r) = r + qoar^ + cr^ + dr'' + er^ + ©(r"^) 

and 

fk+qo = Vk + 9000 + qobrl + 0{rl) 
which we proved in Lemma [5] and 

cos{^k+qo - (^fc) = 1 - ^b^q^rl + 0{rl) 

and get 

zl^{a'+b')qlrl + Oirl). (48) 

It is easy to show that the behaviour of y^. is the same as for fixed and 
rational rotation angle (Lemma 3). The sequence yk is 

yl = rl+ rl^^ - 2rkrk+i cos(v3fe+i - ipk) = 

= rl+rl^^-2rkrk+iC0s{Q^ + hTl + 0{ri)) (49) 

and the behaviour of cosine of the rotation angle is 

cos{Q^ + brl+0{rt)) = cos{Qo) - {brl + 0(4)) smiOo) - 

= cos(eo)-6r2sin(eo) + 0(r4). (50) 

It follows 

yl = rl + {rk + arl + 0{rt)r- 

~ 2rk{rk + arl + O(r|))(cos(eo) - brl sm{Q^) + 0{rt)) = 

= 2rl + 2ari + 0(4) - 2rl cos(eo) + 2brt sin(eo) + 0(4) = 

= 2(l~cos(eo))r^ + 2(a + 6sin(eo))r^ + 0(rf). (51) 

The term Wk is of a form 

^k = ^fc+l + ^fe+go - 2rfe+irfe+,„ cos{tpk+qo - <^fc+l). 
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We use that 

^k+qo = + (<Zo - l)eo + (go - l)brl + 0{rl) 
which we proved in Lemma [5] and 

cos{ipk+qo - "Pk+i) = cos((go - 1)00 + (go - + C(?-fc)) = 

= cos((go - 1)60) - (6(go - l)rl + 0{ri)) sin((go - l)eo) - 

- (5(go-lH + C^(rDf ^°^((^"-^^^°^ 0(.|) = 

= cos((go - 1)60) - fe(go - l)rl sin((go - l)eo) + O(r^). 

It follows 

wl = {ru + arl + 0{Ti)f + {rk + q^arl + 0{ri)f-2{rk + arl + 0{Ti)) 

(r-fc + goar^ + O(r|))(cos((go - 1)60) - (go - l)hrl sin((go - 1)60) + 0{rl 
= 2rl + 2qoart + 0(r|) - 2rl cos((go - 1)60) + 2qobrt sin((go - 1)80) + 0(r 
= 2{l-cos{{qo-l)eo))rl+0{4). 

We put the behaviours of {r^} in the expressions for yk, Zk and Wk, and the 
claim is proven. The last claim will be showed analogously as in the proof of 
Lemma 3, only with different coefficient va? + 6^. I 

Remark 4. For the case with a > 0, we study the inverse system. The 
lemma is also true for odd a > 3. 

By direct calculation we will show in the next theorem the same claim as in 
the Theorem 3, but for the general rational case and a = 3. 

Theorem 4 Let 

r ^ r + ar^ + Oir'^) = f(r) 

LP ^ ip + Oo + br^ + 0{r^) = g{ip, r) (53) 

be a system with a ^ 0, Qq = for some p, g g Z and e^^° ^ a/T for 

n — 1,2,3,4. Then there exists ri small enough such that the discrete spiral 
r(ro, (po), for To < 7'i is Minkowski nondegenerate, and dims r(ro; lySo) = 2/3. 

Proof. 

Without loss of generality, we may assume that a < and p = 1. Namely, 
if a > 0, then xq is unstable fixed point. In jl2j . Lemma 3.1 we showed that 
the inverse map of f{r) = r + ar^ + 0{r'^) for a > is of a form f^^{r) = 
r — ar^ + 0(r*), so we are back to the case o < 0. 

The proof will be done for the case of fixed displacement of angle, that is, 
g((p, r) — (p + Oq. In the previous lemma it was shown that in the general case 
the behaviour of overlapping is the same, and that all other preconditions are 
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valid as well. The difference in the quotient of the behaviour of {zk} can have 
an influence only on the explicitly calculating Minkowski content, what we will 
not calculate here. 

It follows from the previous lemma the the first overlapping is Zk < 2e if 
k > Kq for some fixed Kq. But if we denote the discrete spiral as 

r(ro, fo) = ^^Jo^Tiin, ifii) U T{rKo,^Ko), 

where Tj are the points with polar coordinates {ri,ipi), and T{rKo,'PKo) is a 
discrete spiral with the initial point {rKo^fKo) > that is, 

Because of the finite stability, it holds that 

dimBT{ro,(po) = dimBT{rKo, VKo)- 

Let £o be such that the least number mi{so) for which Zk < 2eo is greater 
then Kq. We take fixed e such that e < eq- Now for such e we try to find the 
least natural number mi (e) such that Zk < 2e, and we get that 

mi(e) ~ {qoa)^s~i, 

that is, 

2 2 

^i£~3 < mi(e) < Bie~3. 

The second overlapping is the overlapping after which we are in the core since 
we see that the overlappings Wk < 2e and yk < 2e are simultaneous, and we get 

m2{e) ~ 2(1 - cos(eo))£"^ 

that is, 

A2£-^ < m2(e) < B2£"^. 

Now we will describe the £-neighbourhood of discrete spiral r{rKo,VKo)- 
As has three parts: 

1. \Ai\^ - £-neighbourhood from initial point until first overlapping (taili) 

2. 1^2!^ - £-neighbourhood Minkowski sausage from first until second overlap- 
ping (tail2) 

3. 1^43]^ - £-neighbourhood from second overlapping until the fixed point (core) 

Now we estimate part by part. So, 

= (mi(£) - Ko)e^TT -(£"§- Ko)e''Tr ~ ei , (54) 

that is, 

Ai7r£3 - Kos'^TT < \Ai\^ < Bme^ - Koe'^TT (55) 
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Now the part from first until second overlapping is 

m2{s) — l 



1^2 L = ("22(e) — wi(£))e^7r — £^ (2 arccos(^) — — sin(arccos(^))) 

^ — ' 2e £ 2£ 



k—mi (e) 

We use the equality 



(56) 



sin(arccos x) = \/l — x-^ 
and the following sequences which are valid for < 1 of a form 

TT" S ^ In A-^ 

arccosx — — — a^x — aix — a2X — ... — a„.T — . 
with ao = 1 and a„ = 2-4'6.!(2i^)"2n+i) for n = 2, 3, . . . and 



Vl-x^ - 1 + bix^ + b2X^ + ... + 6„2;2" + 



1 , l-3-5...(2n-3) 

where h = 2.4.6...(2„) ■ 

So we simplify the sum 

m2{e)—l 

E(2arccos( — ) ^ sin(arccos( — ))) = 

fc— mi (e) 

m2{€)~l 3 2n+l 

- T. p<i-|-°.3=----^)--- (^v 

k=mi(e) ^ ' 

E/ _ „ffc _ «! + Ol 3 _ «2 + O2 5 _ _ a« + On 2n+l n 
4^3 23£5 22"£2»+l^'= 

k—mi (e) 

Now we put it in ([56| and get 

m2(£)— 1 m2{e) — 1 

\A2\^^im2ie)~mi{e))e'^'n:~e'^TT ^ 1 + 2£ ^ Zfe + 

k—nii (e) k—mi {s) 

m2{e)~l m2(£)-l male)-! 

+^ E ^.^ + ^ E 4 + ... + !^ E ^r^ + -(58) 

/c— mi(£) k—mi(e) fc— mi(e) 

Notice that the first two monomials are abbreviated, so by including the expres- 
sion 

an + On = 



2n+ 1 
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we have 



™2(e)-l , m2(e)~l m2{e)-l 

Ai^mi (e) k—nii (e) k—jni (e) 



7712 (e) — 1 

T E + (59) 



(2n + l)(2e)2» 

K=mi (e) 

Now we only need to find a "good" estimations of these sums. 
If Azk^2 < Zk < Bzk^^, then 

6-1 6-1 6-1 



k—a k—a k 



We will approximate the sum of powers of k by using the integrals, since it is 
easy to show that for the decreasing function / it holds 

b-l „6-l 



/ f{x)dx < J2 fik) < I fix)dx, 



and we get 



So, the lower bound for the sum is 

k—mi (e) 

while the upper bound is 

9R2n+l i i i i 

y 2„+l < ^ [(J_)3n+|£2„+i^/3^^i)(J_)3«+|^2„+l_._L)3n+i^6«+l^(5. 

^ ^ - 6n + l Ai ^ 2 Ai 

/c— rni (e) 

By using this bounds for power sums of z^, we estimate (|59p and get 

It is easy to show that < 1 so the sum 

" 7^ + 7^ it (2n + l)(6n + 1) |^ ) ^ ^ 
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converges. But we want to show that 

Ci >0, 

that is, 



bn 



A. 



(2n+l)(67i + l)^2/B3^ A,^' 



> 



Recall that 6„ < 0. We know that for x 



< 1 



(61) 



n— 1 ^ ^ ^ ^ 



= VI -a;2 - 1 > -1. 



ri=l 



Since '^^'^j^ > Ij the inequality (pT|) is true. 



Now it follows that 



where Ci is given by and C2 = 4B2((^)5 - -^). Then 



The upper bound for the nucleus is 



(62) 



(63) 



(64) 



So we get the estimation for the Minkowski sausage in the rational case, and it 
follows fromdSi]), §^ and ([M]) 



-2-s 



that is. 



\Al 



Now we have 
and 



< (BiTT + Ci)£^-3 + (C2 + C3 - ^)e' + 0(£"+3 ) 



7W*3(A) < BiTT + Ci < 00, 



dimsF < -. 



But for the lower bound we have 

I^L |AiL + |yl2L + |A3| 



-2-s 



-2-s 



It holds 



Mi{A) > MtiAi) + Mi{A2) + MfiAs). 



(65) 
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Now for d = I wc have MS {Ai) > Am > 0, MS {A2) > and MS {A3) > 0. 

2 

So we get MS (A) > 0, and 

2 

dim„r > -. 

Finahy, it is proven that F is a Minkowski nondegeneratc and dim ^ F = | . ■ 

Remark 5. It is possible to generalize the claim from the previous proposi- 
tion for a = 2k + l with the change of the nonrcsonant condition into e*®° 7^ \/T 
for ri = 1, 2, . . . , 2fc + 2, and then dim^ F(ro, ipo) = 1 — ^. 



3.4 The unit-time map of Hopf-Takens bifurcation 

Since the order of overlapping in the irrational case is even more complicated, 
the direct calculation of the box dimension will be more demanding then for 
the rational case. In order to avoid it, in this case we are using the connection 
between the Hopf and Ncimark-Sacker bifurcation through the unit-time map. It 
means that we will prove the box dimension result only for the discrete systems 
which are unit-time maps of continuous systems with Hopf bifurcation. 

At the beginning we have to demonstrate how the unit-time map of Hopf- 
Takens bifurcation looks like. We consider the continuous dynamical system 

± = F(x) (66) 

where x G M" . The simplest way to extract the discrete dynamical system from 
(|66[) is by using the flow of a system (j)t{x). Namely, we fix tg > and wc observe 
the system on X generated by the iteration of flow (j)tg (the map with shift 
along the trajectory of ([66]) ). Now we chose 1^ = 1 and get the discrete system 
generated with the unit-time map 

x^(t>i{x). (67) 

It is easy to show that the isolated fixed points of (|67l) corresponds to the isolated 
singularities of (|66|) . It is also known that the corresponding singularities and 
fixed points are simultaneously hyperbolic or nonhyperbolic. 

Now, in order to find the unit-time map for the generalized Hopf bifurcation, 
we start with the normal form for the planar system with the two complex 
conjugated eigenvalues iwi which has a vector form 

-ujy + {x^ + y'^){ax - by) 

u)x + {x^ + y'^){bx + ay) (69) 



X = 

i.e. by components 



-LO 

UJ 
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or in polar form 



r = ar^ 

ip = uj + br"^. (70) 
For the system (|70| . the third Lyapunov coefficient is 

9^ 
= — a. 

So if a 7^ 0, then V3 7^ 0, and by Theorem 5 from 31 it follows that dims F = 
I where F is a spiral trajectory near weak focus. Moreover, F is Minkowski 
nondgenerate, that is, 

< Al^(F) < X*'*(F) < 00 (71) 

with d = |. This result will be used to get the box dimension of Neimark- 
Sacker bifurcation. As we know, the unit-time map can be easily obtained by 
the method of Picard iterations. We get 

<pi(x)=e^x+||x||^(^ ^ -^)e^x + 0(||x|r) (72) 

where 



cos it} — sm w 
sinw cosw 



(73) 



In the complex form it is 



/(z) = e*"z + |z|' (a + ih)e'^z + 0{\zt). (74) 
The polar form for / is 

r ^ r + ar^ + 0(r'^) 

if Lp + uj + br'^ +0{r'^). (75) 

Now we can conclude that the normal form map for the Neimark- S acker 
bifurcation can be approximated until the 3rd degree by the continuous system 
which exhibits the Hopf bifurcation. We can also see that all the unit-time 
maps of the Hopf bifurcation has the form / with the orbits which lie on the 
trajectories of the system (|69p . See Figure 6. By using the property of box 
dimension regarding the subsets (A C _B ^ dims A < dimsB), we get the 
same upper bound for the box dimension as before, that is, |. Now we would 
like to prove that the lower bound is the same, in order to prove that the box 
dimension is |. 
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Figure 6a Figure 6b Figure 6c 

a = -1, b = l, Lu = 1 a = -1, b = I, Lo = ^ a = -1, 6 = 0, w = f 

Now we observe the generalized planar continuous system 
X — —Loy+{x'^+y'^)^{ax~by) 

y ^ ojx + {x'^ + y'^)^ {bx + ay) (76) 
which in polar coordinates is 

r = ar^''+^ 

(fi = uj + br^''. {77) 
Lemma 5 // we have the planar system 

(p = u + br^'', (78) 
then the corresponding unit-time map is of a form 

r ^ r + ar^'^' +0{r^''+^) 

if ^ (fi + uj + br^'' +0{r^''+'^). (79) 

Remark 6. The lemma can be easily proven using the Picar iteration 
method. 

Remark 7. By adding the terms of higher order to the continuous system, 
the form of the unit-time map doesn't change. 

The following lemma is a generalisation of the box dimension result from 
(Theorem 9(b)) and it can be shown analogously. 

Lemma 6 Let T be a spiral trajectory near the origin of planar system 

r = ar^''+^ 

(p = oj + br^''. (80) 
Then T is a Minkowski nondegenerate and we have 

dim. 2(1-^). 
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Proof. 

This lemma wih be proven by using Theorem 7 from |28) . It means that we 
have to show that r = f{f) is a radiaUy decreasing function, and that it satisfies 
the following conditions: r = f{ip) ~ <fi^^ , f{ip) — f{ip + Stt) ~ ip~^~^ and 
~ (p~2k-^ and < M^ip^^k, From the assumption a < 0, it 

follows that r ~ f{if) is a radially decreasing function, that is, the origin is a 
stable focus. 

Now we seek for the solution of a system (|80l) . and get 

f ar^k + 1 



or 



By integrating, it follows 



a{2k) a 

ie. 



' - ln(r) =^ + C, 



a(Zk) a 

It is easy to show that (p = $(r) ~ 7--2fe^ g^j^^^ then we have r = f{(p) ~ ip^^ . 
Then \f'{ip)\ ~ ip~^~'^. By Mean Value Theorem we get f{(f) - f{ip + 2n) ~ 
(p^2fc-i. We know 

so it follows the last condition on the derivative 

\f"{ip)\ = |/'(v3)2$"(r)$'(r)-i| < ap-i-^if-^-^if^+i < A^ip-^ . 
Now we have dims F = 2(1 — jFTt)- ■ 

Until now, we observed the truncated normal form for the Hopf bifurcation at 
the bifurcation value — Q because then the origin is a nonhyperbolic singularity 
and the box dimension is nontrivial. But, if we consider the whole one-parameter 
family of planar vector fields which exhibits the Hopf bifurcation aX ^ — 0, it is 
of a form 

X = d^x — [u + cii)y + {ax — by){x^ + y^) 

y — [uj + c^)x + djiy + {bx + ay){x'^ + y"^) (81) 



or in polar coordinates 



r ~ d^r + ar^ 

(f = w + c/i + (82) 
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with the transversality condition d ^ and nondegeneracy condition a 7^ 0. 
Analogously as before, we get the unit-time map 

r I— > r + d^ir + ar^ + 0{^^r, ^r^ ,r^) 

ip + Lo + cfi + br^ +0{n^,^lr'^,r'^). (83) 

It is known that the transversality condition for the Neimark-Sacker bifurca- 
tion is also d 0, and the nondegeneracy condition is also a ^ 0. Since the 
nonresonant condition for the Neimark-Sacker bifurcation ensures that the nor- 
mal form is given by (|83l) . we don't need this condition since every unit-time 
map is of a given form. So, in the same time the system (|8ip experiences the 
Hopf bifurcation for fi = and the corresponding unit-time map experiences 
the Neimark-Sacker bifurcation for ^ = 0. At Figure 7 we can see the Hopf bi- 
furcation and Neimark-Sacker bifurcation with lu = 1, a = d = l and b = c — 0. 




Figure 7a fi > Figure 7b fj, — Figure 7c /i > 



3.5 Irrational case 

In this chapter we will calculate the box dimension of the orbit of discrete system 
which exhibits the Neimark-Sacker bifurcation in irrational case. The normal 
form at the bifurcation value is 

(p ^ ip + Qo + br^ + 0{r^), (84) 

where 9o = 27r/3, (3 G K\Q (irrational case). 

Before we begin with the ordering of overlappings, notice that the angle 
mapping depends on the radius. So we observe the expressions for yk, Zk and 
Wk- Since is irrational, then there exists the lowest natural number such 
that go9o > 27r, that is go = ^ ■ Let us look at the behaviour of the sequences 
{ijk}, {zk\ i {wk\ defined in Subsection 3.2. 

Lemma 7 Let 

rk+i=rk + arl + 0{rl) 
(pfc+i = (^fc 60 + brl + 0{rt) 
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be a given system with a ^ 0. Then for Oq = 27r/3, /3 G M\Q where qo — ^ 
we have 

Vl = rl + rl+,-2rkrk+icos{Qo) + 0{rt) = 2{l-cos{Qo))rl+0{rt) 
4 = rl + rl+^^-2rkrk+,,cosiqoeo) + 0{rt) = 2{l-cos{qoeo))rl+Oirt) 
wl = rl^^+rl^^^-2rk+irk+q„cos{{qo-l)Qo) + 0{rl) = 
= 2(l-cos(go-l)eo)r2+0(r4). 

Proof. 

The behaviuor of {yk} can be showed as in Lemma 3. We consider 

zl = rl + rl^q^ - 2rkrk+qo cos{ipk+qo - Vfc)- 

Using the expressions for rk+q^ and fk+qo , we get 

4 = rl + {rk+qoarl+0{rl)f- 

- 2rk{rk + qoarl + O(r^))(cos(goeo) - <Zo&r^. sin((7oeo) + 0(r|)) = 

= 2r2 + 2q^art + 0(r^.) - 2rl cos{qoQo) + 2go?'r^ sin(goeo) + 0{rl) = 

= 2(l-cos(goeo))r^+2qo(a + 6sin(goeo))r^+0(r|). (85) 

Analogously we get 

= {rk+arl + 0{ri)f + {rk+q^arl+0{rt)f-2{rk + arl + 0{ri)) 

(rk + qoarl + O(r|))(cos((go - 1)60) - {qo - l)brl sin((qo - 1)60) + 0(r 
= 2rl + 2qoart + 0{rl) ~ 2rl cos(((7o - 1)60) + 
+ 2((jo-l)6r4sin((go-l)eo)+0(rf) = 
= 2(l-cos((Qo-l)eo))r2+0(r4). 



Before we start to calculate the e-neighbourhood of the orbit of the 
system (j84p , we are interested in the order of overlappings of neighboring terms 
of sequences y^, Wk and Zk since all of them behave like fc~ 2 but with different 
coefficients. If the first overlapping is between the adjoining points (yk < e), 
then the lower bound for of that overlapping is in fact the £*-neighbourhood 
of the corresponding spiral trajectory with the radius e* > ^^^ e. That will 
be the part of the tail for the £*-neighbourhood |^e*| for spiral trajectory. 
Furthermore, if for the spiral trajectory with e* the entering in the nucleus is 
after the overlapping yk < s, then that nucleus with e* is the lower bound for 
the nucleus with e. 

Therefore, if the first overlapping is the overlapping of adjoining points, then 
we have the result as we will see in the proof of the theorem . Now we want to 
show that this is always true. 
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Lemma 8 Let be given a system with Qq = 2TTf3, p G Then there 

exist 01 small enough and Kq big enough such that for all Qq < &i and k > Kq, 
it holds yk < Zk- 

Proof. 

By considering the behaviour of the sequence 

yk = d{Ak,Ak+i) = \/rl + r^+i - 2rkrk+i cos{Qo) + 0{rl), 
we see that 

yl = (rk - rk+i)^ + 2rkrk+iil - cos(eo)) + O(r^) > (r^ - rk+i)^ + ^. 

Notice that if the angle 9o is smaller, then yk is closer to (r^ — f/c+i) + fr, ie. 
when 00 — > 0, then yk — > {rk — rk+i). By considering the behaviour of 

Zk = d{Ak,Ak+q„) = ^rl+ rl^g^ - 2rkrk+q„ cos{qoQo) + 0{rl) 
we notice that 

4 = in ~ rk+q„ f + 2rkrk+q^il - cos(goeo)) + 0{rt) > {rk - rfc+gj^ + ^. 

We see that if Qq is smaller, then (7o0o(niod27r) is smaller, so Zk is closer to 
Tfc — rk+qo + ji- But also when Go — > 0, then go ^ c» and — rk+q^ increases 
from rk — rk+i + p- to r/j + p-. It means that we get: if Qq 0, then yk 
decreases to rk ~ rk+i + p-, and Zk increases from — rk+i + p- to + ||. 
Hence, we can chose 9o small enough such that yk < Zk, for every k big enough. 
Let 6i be such an angle, ie. for all 9o < 9i we have yk < Zk, for k big enough. ■ 

Theorem 2 showed that 

2 . 4 
— < dimoF < dimnr < — . 
3 ^ - 3 

It is already proven that for 9o = 27r/3, /3 G Q the box dimension is equal to the 
lower bound. In the next theorem we will prove that in the irrational case the 
box dimension is equal to the upper bound. 

Theorem 5 Box dimension for Neimark-S acker bifurcation - irra- 
tional case 

Let we have a discrete planar system 

r r + ar^ + 0{r^) = f{r) 

(f ^ f + Qo + br^ + 0{r^) = g{r,ip). (87) 
which is the unit-time map of the continuous planar system 

r = ar^ + 0{r^) 

ip = eo + 6r2 + 0(r4), (88) 
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and let = 27r/3, with (3 irrational. Then there exists ro > such that the 
sequence S{ri,(pi) = {{rk,(pk)} defined by rk+i = f{rk) and tpk+i = g{rk,ipk) 
where ri < r^, is a Minkowski nondegenerate , and we have 

4 

dmiB S{ri,^i) = -. 

Proof. 

First we should check the ordering of overlappings. If w < Qi, then it 
follows from Lemma 8 that yk < for k big enough. So in that case the 
first overlapping is for the sequence yk- But, what if the condition uj > Qi 
isn't satisfied. Then we would like to find some other system which is also the 
Tg-time map of (|88p . but with less then Oi. And, of course, we would like 
that system has the same box dimension as the initial system. We consider the 
system 

r ^ fi (r) 

V ^ 9i{r,ip) (89) 

where /i and gi are such that ff = f ^ gf = g for the lowest natural number q 
such that ^ < 6i. It is easy to see that the previous system is of a form 

r h-^ r+-r^+0{r^) = ,f{r) 

q 

if ^ <^+^ + V+0(r4) =g(r,(^). (90) 

q q 

Furthermore, notice that since the system (|87p is the unit-time map of the 
system ((55)) , then in fact the system (lUHl) is a i-time map of the same continuous 
system. It remains to show that the box dimensions of the orbits of systems 
((87|) and ([90)1 are the same. Namely, if F* is a discrete spiral of ((90)) with the 
initial point (ro,i^o), then 

F*(ro,^o)=ULiF,;(r^,^^), 

where F^ are the discrete spirals of (|87)) with the initial point 

irh,^l)^{fr\ro),gl-\ro,^o))- 
Now, it follows from ([IJ 

dims F*(ro, (po) = dims F(ro, (po), 
where F is a discrete spiral of ()87|) with the initial point (ro, (po)- 

Hence, without loss of generality we may assume that cj < Oi, that is, that 
the first overlapping is yk < 2^. Furthermore, since we would like to find the 
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lower bound, we can take the inequality yk < £ instead of yk < 2e. The difference 
in the estimation of e-neighbourhood will not influence the lower bound. Then 
we get the lower bound for the corresponding spiral trajectory by using the 
"annular" spiral. See Figure 8. If yk < e, then the radius of "annulus" is e* 
such that £* > ^'^^e. 

We denote hy A = r(ro, (po) discrete spiral, and with T spiral trajectory of the 
system (|88l) . 



Figure 8 

We divide e-neighbourhood into several pieces to simplify the estimation 
process. 

PART 1. j4i.£ = £-neighbourhood from (ro,<y5o) until the first overlapping 

= £2^(mi(e)-l)>£27r(cid£-i-l) = 7r(cid£^-e^) (91) 
^ > Acid£'-^-e') (92) 

where mi{e) > cid£~i since yk > [vk - rk+i) ^ 

PART 2. A2,e ~ £-neighbourhood from the first overlapping until the end 
consist of two parts: the tail |A2,e| which contains all possible overlappings and 
the nucleus l^s.el- This part will be estimated by the £*-neighbourhood with e* 
from the corresponding spiral trajectory of continuous systems near focus which 
also contains two different parts: the tail Ir^^e*] from </'mi(e) do ip2{e*) (this is 
the angle after which we are in the nucleus) and the nucleus |rj £_,|. Hence, it 
holds 

\A2.e\ + \A^,e\ > \Tr.,e'\ + \T,,,J (93) 



Notice that the right side of £-neighbourhood is continuous spiral F for g 

4 
3 



[(/Jmi(6)j oo). Now we use that F is a Minkowski nondegenerate and for d — 



we have 

X^(F) = liminf^>0. 

e-^0 £2-d 
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It follows 



Now we include s = d to get 

""^'>™33f^- (95) 



l^2,el + l^3,el ^ |r.,.. | + |r,, , 



It follows from e* > that e'^ ^ > ( ^^ ") (e*)'' ^, so we get 



l^2,el + l^3,sl ^ ,y3-i,2_jr,,,.| + |r,-. 



^(^)^'^ ' VO^V- - (96) 
Now, taking into consideration the parts 1. and 2. we have 

MtiA) ^ hmmf + l^^/:-' + > 0+ (^^) ^ 'A^S?(r)>0. (97) 

Therefore, we have the lower inequality for the e-neighbourhood. 

On the other hand, it is quite obvious that the e-neighbourhood for the 
discrete spiral is contained in the e-neighbourhood of the corresponding 
continuous spiral Tc^e ie- 

+ \A2,e\ + \A3,e\ < \T,t,e\ + l^cnj = \Tc,e\ (98) 

where Tct,e is a tail of e-neighbourhood of continuous spiral F, and Tcn,e is a 
corresponding nucleus. So we have 



M*'{A) = hmSUp + +1^221, ^ ^ (gg^ 

0<X^(A) <7V/(*'*(A) <c5o 



Hence, we have 



and it follows that 

dims (A) = d. ■ 

Remark 6. Using the Lemma 5 and Lemma 6, it is possible to prove 
the analogous result for the generalized Neimark-Sacker bifurcation, that is, 
Chenciner bifurcation (see [S], [5], P^,)- If we apply the obtained result for the 
box dimension of the Neimark-Sacker bifurcation in irrational case to Chenciner 
bifurcation, we get that the box dimension of the orbit near the origin at the 
bifurcation point (/3i,/32) = (0,0) is dims S{ri,ipi) = 2(1 — |) = |, as we can 
see at Figure 9. 
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Figure 9a f{r) = r — 
g{(p,r) = 1 + 



Figure 9b f{r) = r — 
g{ip, r) — 0.5 + 



Remark 9. Regarding the Neimark-Sacker bifurcation, we can conclude 
that the difference in the box dimension is connected to the different order of 
overlappings. Namely, in the rational case it is Zk < yu for all k big enough, 
while in the irrational case the ordering is reverse. We can also see that the 
Neimark-Sacker bifurcation has a nondegeneracy condition of order three. So, 
like in all previously studied bifurcation (see [12],[28]), the box dimension is 
connected with the order of nondegeneracy. But in the case of this bifurcation, 
it is also connected to the rationality of displacement angle. Moreover, notice 
that the discrete spiral in rational case has the same box dimension as Poincare 
map of continuous spiral near weak focus (see |31p. while in the irrational case 
we get the box dimension of continuous spiral. It is also interesting to notice 
that the box dimension of bifurcated object differs from (in the case of periodic 
orbit on the invariant circle) to 1 (in the case of dense orbit on the invariant 
circle) . 

Remark 10. It will be interesting to further explore Neimark-Sacker bifur- 
cation of limit cycle and the bifurcation of the systems in higher dimension with 
Neimark-Sacker bifurcation such as fold-Neimark-Sacker (fold-NS, see [3] ), flip- 
Neimark-Sacker (flip-NS, see [16]) and double Neimark-Sacker (NS-NS, see [16]). 

Remark 11. This connection via unit-time map between the discrete and 
continuous dynamical systems can be applied to fractal analysis of nilpotent 
singularities of planar vector fields. See (T4] . 
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